Effective action and motion of a cosmic string 



Malcolm Andersontl 

Department of Mathematics, Edith Cowan University 
2 Bradford Street, Mount Lawley, WA 6050, Australia 



- 1 — \ 

X 



Filipe BonjourQ and Ruth Gregory 
Centre for Particle Theory 
University of Durham, Durham, DHl 3LE, U.K. 

John Stewart! 
Department of Applied Mathematics and Theoretical Physics 
Silver Street, Cambridge, CB3 9EW, U.K. 
(February 1, 2008) 

I> 

' We examine the leading order corrections to the Nambu effective action for the motion of a 

, cosmic string, which appear at fourth order in the ratio of the width to radius of curvature of the 

string. We determine the numerical coefficients of these extrinsic curvature corrections, and derive 
the equations of motion of the worldsheet. Using these equations, we calculate the corrections to the 
motion of a collapsing loop, a travelling wave, and a helical breather. From the numerical coefficients 
we have calculated, we discuss whether the string motion can be labelled as 'rigid' or 'antirigid,' 
and hence whether cusp or kink formation might be suppressed or enhanced. 
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■ I. INTRODUCTION 

cn 
i> . 

, The study of topological or vacuum defects is of importance in many areas of contemporary physics. In high 
1 energy physics, a defect will generically occur during a symmetry breaking process where different parts of a medium 
On ' choose different vacuum energy configurations, and the non-compatibility of these different vacua forces a sheet, line, 
(-H : or point of energy where these non-compatible vacua meet. The relevant vacuum order parameter then becomes 
O indeterminate — this is the defect. A defect may be topological ||l[], in that it is the topology of the vacuum that 
' simultaneously allows formation, and prevents dissipation, of these objects — but other types of defect are also possible. 

For instance, a defect may be stable dynamicalhv' (i.e. classically, due to energy considerations) but not topologically, 
(— I ' as it happens for semilocal |^ or electroweak [g| defects. A defect can even be 'topological' and unstable, as in the 
case of textures [Q, but nonetheless of physical importance. 

In cosmology, there has been much speculation that topological defects (stable or otherwise) might have played an 
important role in structure formation Q . In general, there are two main concerns when considering the cosmological 
$^ ' effects of topological defects: their gravitational effects and their dynamics. Any theory concerning galaxy formation 
must be able to allow or constrain the presence of strongly self-gravitating objects. But the dynamics of the defects 
are in fact even more important, for it is the dynamics that determine to a large extent the shapes of the gravitating 
defects. For instance, if cosmic strings did not intercommute, any network would rapidly become tangled and would 
not obey a scaling law; such a configuration would be in conflict with the universe we see around us today. Even with 
intercommutation Q, strings which are strongly rigid and hence straight will have different gravitational effects to 
those that are very crinkly ||^,^ . 

The dynamical behaviour of a defect is generally assumed to be approximated by an effective action, a description 
which models the rather large numbers of degrees of freedom of the full field theory by the smaller number of degrees of 
freedom based on the position of the core of the defect. Attempts to derive effective actions or equations of motion for 
topological defects have commonly focussed on the strong coupling limit, meaning that of large values of the coupling 
coefficient A of the relevant Higgs field. In this limit, the defect becomes infinitesimally thin and effectively decouples 
from the other (infinitely massive) particles in the field theory. The study of the effective motion of topological defects 
has been extended away from the limit A ^ oo to cases for which the thickness is small but not exactly zero. 

The resultant effective action generically contains a 'zero-thickness' term proportional to the area of the defect p2| , p^ , 
and extrinsic curvature terms which appear at quadratic order in the thickness |^-|llj • The exact form of these second 
order terms is dependent on whether an 'off-shell' or 'on-shell' approach has been used to derive the effective action 
as was explained using the example of the domain wall in p[ , and one finds that in a self-consistent order by order 
solution of the equations of motion, the only quadratic correction appearing is due to the geometry of the defect 
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worldsurface, and is proportional to its intrinsic Ricci curvature [Q. For the domain wall, this term gives corrections 
to the motion, however for the string this term is a topological invariant — proportional to the Euler character of the 
worldsheet — and hence does not give any correction to the Nambu equations of motion. 

The leading order corrections to the motion of a cosmic string were shown by one of us, |]l5| , to appear at quartic 
order in the string width. In this paper, a systematic expansion of the geometry and field equations clarified an earlier 
discrepancy concerning second order 'twist' terms and derived the fourth order action of the string, which is 
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where 7?,, JC^^u are the Ricci and extrinsic curvatures of the worldsheet (to be defined in the next section), and the 



are numerical coefficients depending on the field theory modelling the vortex (to be defined in section III for the 



Abelian-Higgs model). Unlike the domain wall, whose background solution, and extrinsic curvature corrections can 
be expressed analytically |l^, the prototypical cosmic string solution, the Nielsen-Olesen vortex ||l^, does not have 
a closed analytic form, hence these coefficients must be evaluated numerically. 

As soon as one calculates corrections to the Nambu action, it becomes of general concern whether or not these 
corrections cause the motion of the defect to be 'rigid' or 'antirigid.' The implications of extrinsic curvature terms 
for string motion have been well-explored [l0|Jl^ -[l9| , however, it is not always easy to decide a priori (especially for 
the fourth order terms) whether the strings will be rigid or not, or indeed what one means by rigid p9[ |. 

In this paper we address these issues. We determine the numerical values of the Ui for the Abelian-Higgs model by 
solving the perturbed field equations for the non-flat worldsheet. We then derive the fourth order equations of motion 
from Anderson's action and calculate corrections to three 'test case' trajectories: the circular loop, the travelling 
wave, and the helical breather. These three solutions display different characteristics which should be mirrored in the 
corrections to the Nambu motion if the rigidity of the string is to be determined. The loop collapses to a point [ pO[ , 
and rigidity would be indicated by a retardation of this collapse. A travelling wave on the other hand has been shown 
to be an exact solution of the full field theory |^ , and should not exhibit any corrections if our approach is to be 
trusted. The helical breather (see e.g. |^) is a time dependent solution which is never singular. Rigidity for this 
trajectory is more subtle, since the helical breather is never singular, however, we could call the string rigid if the 
tendency of the correction is to lower the magnitude of the scalar curvature of the worldsheet. As we will see, it is 
rather difficult to give an intuitive criterion for rigidity. 

The layout of the paper is as follows. In the next section we review the formalism for the derivation of the effective 
action. In section III we rederive Anderson's action, and present new numerical results evaluating the coefficients 
appearing in the action. In section |^ we derive the equations of motion of the fourth order string, and in section ^ 
we calculate corrections to three test trajectories. In section VI we discuss the question of rigidity and conclude. 



II. DERIVING THE EFFECTIVE ACTION 

In this section we review the formalism required for the derivation of the string effective action. This formalism is 
largely based on the contents of (l^Jl^ • The problem of building the effective action is to reduce some four-dimensional 
field theoretic action integral 



to some two-dimensional worldsheet integral 
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We follow the canonical approach by setting up a coordinate system based on the worldsheet, and expanding the 
action and equations of motion around this worldsheet. By 'expand' we mean that we do not expect in one go to be 
able to solve the full equations of motion and integrate out (otherwise why find an effective action) but that we will 
be able to express the equations of motion for the system in terms of an expansion in r^K, where rs is the radius of 
the string, and k, a typical scale of the extrinsic curvature of the worldsheet. We then solve the equations of motion in 
the directions perpendicular to the worldsheet, and integrate out the four-dimensional action (^ over these degrees 
of freedom to get an action of the form (^. 

We will look at the particular case of a U(l) local string in flat spacetime (signature -f , — , — , — ). This is a vortex 
solution of the Abelian-Higgs model 



2 



1 ~, 



A, 



(4) 



where V,. 



ieAfj^ is the usual gauge covariant derivative, and F^^ the field strength associated with A^. We 



rewrite the field content of this model in the usual way: 



M^n = - 



so that the vacuum is characterised by X = 1. In terms of these new variables, the Lagrangian becomes 
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(5a) 
(5b) 
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where P^j, is now the field strength of P^. By making a compensating gauge transformation in A^j^, we have absorbed 
the unphysical gauge degree of freedom of the theory. For future reference, here are the equations of motion for the 
new fields: 
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-'g^^'X = P^P^X - ^X(X2 - 1) 
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(7a) 
(7b) 



These equations are known to admit vortex solutions, the best known being that due to Nielsen and Olesen — a 
solution corresponding to an infinite straight static string. A vortex solution is characterised by the fact that the 
scalar field phase has a non-zero winding number around an axis corresponding to the zeroes of the Higgs field (X=0) . 
If we choose the string axis to be aligned with the z-axis, then a gauge can be chosen in which the vortex solution 
takes the form 
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in cylindrical polar coordinates, where X^o and Pno satisfy the (generically) second order coupled differential equa- 
tions: 
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and we have introduced the Bogomol'nyi parameter /? — 2e^/A. /3 — > is the global string, and /? = 1 is the critical 
(supersymmetric) string for which the equations of motion factor to two coupled first order equations: 



rX' = XP, 
P' = lr{X'-l) 



(10a) 
(10b) 



Intuitively we expect the solution for a general worldsheet to be very close to the Nielsen-Olesen solution in suitable 
coordinates, and that we might be able to expand the fields, and therefore the action, around this approximate solution. 
Let us now try to make this idea more concrete. 

First we choose a worldsheet-based coordinate system. (Note this will only be valid within the radii of curvature 
of the worldsheet.) We start by coordinatising the worldsheet by a and r as usual. Since we are in flat space, at 
each point on the worldsheet we have a well defined orthogonal flat two-plane, and such orthogonal planes do not 
intersect within the radii of curvature of the worldsheet. Thus we can thicken the worldsheet to a worldblanket by 
deflning cr and r to be constant on such orthogonal planes. Since the worldsheet W has codimension 2, it has two 
associated families of unit normals {n^}. We choose a cartesian parametrisation of the orthogonal planes such that 

(^) ~ ■ Thus {r, (T, define a coordinate system in the vicinity of the worldsheet. In general it will not be 

possible to choose this system to be globally orthogonal, however, we will assume that the normal fields have been 
chosen so that the departure from orthogonality is minimal, i.e. of order of the extrinsic curvature of the worldsheet. 
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Because we no longer have a flat coordinate system, the connection will no longer be trivial, although the curvature 
components should still vanish. In order to examine the form of these relations, we use a Gauss-Codazzi approach |2^ ] 
(see also Q for more of a physicist's approach). For shorthand we denote by {o"^} = 0, 1) the coordinates {r, a}, 
the worldsheet by W and the Minkowski spacetime by Ai. 

The first fundamental form h^^ of W is defined as 

hf^i, ^ g^i, + n^^_,n^^. (11) 

The tensor h^i, acts as a projector onto W, but is still a tensor in the 4-dimcnsional spacetime. For the metric intrinsic 
to the worldsheet, one uses the familiar 

where X^^{a^) are the spacetime coordinates of the submanifold W. This gives an interpretation of W as both a 
worldsheet in spacetime and a two-dimensional manifold with its own intrinsic geometry. One can define the intrinsic 
curvature of W in the standard way, and the extrinsic curvatures, or second fundamental forms, by 

^.M- = hlK^Vpn^a. (13) 

These extrinsic curvatures /C.^j/ measure how W curves in M. Since the codimcnsion of the worldsheet is greater 
than one, we also have a non-trivial normal fundamental form 



where £y is the alternating symbol on two indices. This represents how near to orthogonality the worldblanket 
coordinate system is; it measures the rotation of the in their own planes as one moves around the worldsheet. Note 
that, despite lying tangent to W, o)^ is a gauge dependent object, depending on the choice of the normal fields — an 
S0(2) gauge group. In fact, a)^ is the connection on the normal bundle of W. 
Now, we may write 

= /C;^^ - e^jCui^^^n.^,), (15) 

where the symmctrization is understood to be only acting on indices of the same type, so that the 'j' above does not 
participate in any symmctrization. This now indicates an alternative definition of the second fundamental forms: 

^it^v = ]^^i9tJ^u + eijW(^n^^) = ^Cihf_,^. (16) 

where Ci denotes the Lie derivative with respect to n^. We may now use the Riemann identity in flat space to derive 
the following relation for the extrinsic curvatures: 

Therefore, contracting with g^'^ gives 

£,/C, = -JC^^ylC^" (18) 



In addition, note that 



and 



CiU^iy ^ n.i^.Viyn'^ ^ Cjyeij. (20) 

The system of equations (p"6|-p0|) give the 'equations of motion' for the geometry of the system, which together 
with (^ form the full equations of motion for the string. 

In order to extract a low energy effective action, we will need to expand these quantities and equations of motion 
in terms of the thickness of the string, r^. To do this systematically we define the dimensionless parameter e by 
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where represents a typical radius of curvature of the worldsheet. We then define the zero thickness limit to mean 
e — > with jji (X rf , the energy per unit length of the string, fixed. Note that this may not necessarily mean that the 
string width is zero, since instead the extrinsic curvature could be zero (a flat worldsheet). However, if k ^ 0, this 
limit corresponds to the conventional A — s- cx) limit. 

We now redefine coordinates so that the worldsheet has thickness and curvature of order unity by setting 



and redefining for consistency the following variables 













Pf^i — P^ • 


(Hi) 


become 








= 2e[i^.^i/ - eij 


C K 








CjUJfj, 


= 0, 



and the field equations become: 
1 



P^P^g'^'^X - -X{X^ - 1) 



g^'g'^Fki + {g'^g^' - g''g'^){edAPi ~ OiPa) + eg^^g'^'pAi 



edA 



g^'g''Fki + {g'^'g'^ - g''"g''){dkPB - edsPk) + ey'^^g^^Fsc 



-PX\g'^h+g^^PA) 



g^'^g^'A; + - g^' g^''){edBPi ~ BiPb) + eg'^5^^i^sc 

g^'^g^'Fu + {g^^g^' - g^' g^^){edcPi - diPc) + eg^^g^^Pco 



-l3X\g=^P,+g^^pB). 



(22a) 
(22b) 

(23a) 
(23b) 
(23c) 

(24a) 
(24b) 
(24c) 
(24d) 
(24e) 



(25) 



(26) 



(27) 



Here, Fij is the tensor Fij defined in terms of the rescaled variables. Finally, we also need the explicit expansion of 
the Lagrangian 

C = Xrf I d,Xd,Xg'i + 2ed,XdAXg'^ + e'dAXdBXg^'' + X^P^P.g^^ - 



2/3^^'' 



Fkig'^'g'' + 4(9fePA - edAPk)g"'g''' + 2eFABg 



ik ^jA 



2{d^pA - edApi){d,pB - edBpjW'g^'' - g'^'g'^) 



eFAc 



AB CD 



A{d,PB - edBP)g'''g'''' + eFBog^'^g 



-^{X'-lf 



(28) 
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This now allows us to expand rigorously in powers of e. Note that so far we have made no assumptions about any 
of the fields or their dependence on the coordinates. We have simply rewritten the equations of motion, scaling 
with respect to the physical dimensionful quantities in the problem, leaving the equations of motion in terms of the 
dimensionless parameter e. 

The procedure is as follows. We first consider the e ^ limit, with the metric and fundamental forms taking 
their geometrically defined values on the worldsheet. Since e ^ is also the flat worldsheet limit, we expect that 
the fields will take the Nielsen-Olesen form, which indeed turns out to be the case. In the next section we will go to 
higher orders in e, deriving the corrections to the geometry and fields, and hence the corresponding effective action. 
However, for the moment we conclude this section by proving that the Nielsen-Olesen solution is in fact the leading 
order behaviour for the vortex fields around the worldsheet. 

To zeroth order 

hence all geometrical quantities take their worldsheet values, K.^i^{s^,x^) — K .^i,{s^ ,Q) etc. Hence 



_(lAB 



(30) 



and thus the equations of motion for and P^^ are 



- dAX + (P^ - Pl)X + \x{X^ - 1) = (31a) 

-dAo + PX'^Pj = (31b) 
-d,d,PA + fiX^pA - 0. (31c) 

These are solved by the Nielsen-Olesen solution, which is plotted in figure |l] against r for 13 = 1: 

X„ = XNo(r), P,j = P^ome, P,A = 0. (32) 

Substitution of this solution back into (^) yields the familiar Nambu action. It is the corrections to this action we 
are interested in. 

III. THE CORRECTIONS TO THE NAMBU ACTION 

We now derive the leading order corrections to the Nambu action. In order to calculate quantities away from the 
worldsheet, such as the metric, we must perform a Taylor expansion off the worldsheet 

Q^Q\o + f AQIo + ^eeC^C,Q\o + ^ee^^C.C.CkQlo + ■■■ (33) 

However, in calculating the Lie derivatives of the metric and its determinant some fortuitous cancellations occur. 
First note that 

CkCjK^^, = eCk {K^,^^K^,)pg^P) = (34) 

and also that 

K, K^.K'^^ + 2K^,K;^K^J =0, (35) 

which follows from a trace identity for the 2x2 matrix K |l5| . Hence to all orders, the metric and the volume jacobian 
are given by 

5/^1/ = gtj.u\o + Sex* {K.^^lo - £ijW(^n^^,)|o) + e'^x'x^ (^K^^^^\oKJ^)\o ~ (36a) 
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(36b) 
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This gives all the information on the geometrical contributions to the action. Now let us turn to the field theoretic 
contributions. 

First note that since 
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5S_ 



5X 



9 a I' 



5S 



5P,. 
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(37) 



vanishes by the equations of motion, first order field corrections contribute at second order in the action, and second 
order corrections at fourth order. 

To first order we may read off the equations for the corrections to the fields as 



I) 







-d.F,,, - K^F„,, + pX^P,, + 2pX,X,P,, = 
-d,{ujAx'e,kF^k) - dAP.A + PX^P^A + pxlujAx'SjkPok = 0. 



(38a) 

(38b) 
(38c) 



The first two equations do not apparently lend themselves to a straightforward solution, however, we note that the 
equation of motion for the worldsheet obtained by varying the Nambu action is 



(39) 



which would indicate that the 'driving' terms in (38a, 38b| ) vanish and hence the appropriate solution is X-^ — P^j — 0. 
Indeed, ||T^ demonstrated that unless the trace of the extrinsic curvature vanished on the defect, the first order 
perturbation equations have no solution that is regular and bounded. Thus in fact the appropriate solutions to the 
first two equations of (|3^) are Xi ~ and P-^j — 0. 
The last equation of (|3§|) is solved by 



(40) 



the presence of this term, as pointed out in |I5|, guarantees the gauge invariance of the effective action to worldsheet 
S0(2) gauge transformations, for on substituting in the form of the fields, this correction to Pa exactly cancels the 
'twist' terms of Maeda and Turok, and we arrive at the second order result 



S 



I 



where 



ai = 



27r77' 

2 



rdr X 



r^dri X 



-1 



2"! 



AB 



5/2 



r2/3 

pl2 



r2/3 



1)^ 



(41a) 
(41b) 

(42a) 
(42b) 



For the fourth order action we need the second order corrections to the fields, the volume factor already being exact 
to all orders. From (p5[-p7|) we may read off the equations for these as 



1 



dAX^ + x^K^ABK,d,X^, + P„^P,^X^ + 2F„,P,,X„ + -X^{2,X^^ - 1) = 
-a,F,y + x'^K^abK^^'F,,, + PXIP,, + 2I3X,X,P,, = 

-dtiujAX^EjkFik) - didiP^A 0. 



(43a) 

(43b) 
(43c) 



To simplify these, and to remove the explicit X-dependence we will decompose in cylindrical harmonics by setting 



X. = IgK^AB 



gx'^K^ABK 



AB 



P.<^ = ^I^K^AB + q^X^'K^ABKf'' 
P,r=qreMK,ABKf^x'^, 



(44a) 

(44b) 
(44c) 
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where 



■1^3 



XX- 



This gives two sets of coupled second order differential equations 



P^g , 2P,X,q^ , 1 



.9(3X2 _ 1) ^ 



= -rP' 



and 



= -rX' 



There are no additional corrections to Pa at this level. 

In the critical case /? = 1, the above equations reduce to two sets of first order ODE's: 



rq 



Xo9 + Po 



and 



rP' 

rg' = X„q^ + P,g+—^, 



2g = rqrX„ - 
2qL-4qr = 2rgX„. 



rP' 

0_ 

x„ ' 



(45) 

(46a) 
(46b) 

(47a) 

(47b) 
(47c) 

(48a) 
(48b) 



(49a) 

(49b) 

(49c) 



Note that just as ( [47D is really two coupled second order ODE's, ( [49D is really two coupled first order ODE's. The fact 
that (^, ^) reduce to such a first order form for [3 ~ lis reassuring, since we do not expect to destroy supersymmetry 
simply by bending the string! 

Substituting ([44|) into the action, and after some algebra, one arrives at Anderson's result: 



S ■ 



/ d^a^l 1 



AB 



Q!2 



{J<,ab) + — 



M J J 



where 



Q!2 



as 



■K-q 



dr 



dr 



pi 4^p2 

r^X[^{2g~g) + ^{2q^-q^-rqr)+ " 



/3 



2rX5 + ^(g> + rg,) 



4rP2 

0_ 



(50a) 
(50b) 

(51a) 
(51b) 



If /? = 1, it can be seen from (^, ^) that 2a2 + a^^ —Q. 

To determine a2 and as, we need to solve equations ( |4^ ) and (|47|). This was done numerically by relaxation 
methods, using NAG routine D02GAF |2^]. Sample solutions for both sets of equations for the critical case f3 = 1 
are plotted in figures ^ and |^. The first set of equations can be readily implemented and was solved b y req uiring that 
g{0) = g{oo) = q^{0) ~ 90(00) = 0. Before solving the second set of equations, we expressed Qr from (47c) as: 
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4 + (3r^X§ 



;r ...^2 - (52) 



and eliminated it from the equations, which were then solved by asking that 17(0) — g{oo) = q'^{0) — 50(00) = 0. The 
solutions obtained clearly behave at the origin as: 

9,g,qr^r, (53a) 
Qc/, ^ r'^, (53b) 
^ const, (53c) 

which can be checked to satisfy the asymptotic counterparts of (^ ^7|) . 

We have computed the coefRcients ^, ai,a2 and for a range of values of /3; the results are shown in figures ^ 
and in table |[ The key points to note are that for the range of f3 for which we computed the coefficients (0 < /3 < 100) 
/i,Q;i > and 03, a2 + 0:3 < 0; however, 2a2 + 013 is positive (negative) for /3 < 1 (/3 > 1) and a2 is positive for 
/3 < /3crit — 3.03 and negative for (3 > Pcrit- 



IV. THE EQUATIONS OF MOTION FOR THE STRING 

In order to derive the equations of motion for the string we must express the action ( |50| ) in terms of the worldsheet 
coordinates, with respect to which we are varying the action. To do this, we note that 

K^AB = n^^,A Xi',B = -n^^X^',AB = -n^t.X'^AB- (54) 

Hence, defining 

N^^ = X^-^^,cX^,d'^ (55) 

we see that 

KUb = -X^-'^^'X^.AB = -VVsA (56a) 

K^ABKfK^CDKf'' = X^^'^^X^,cdX''-^'^°X,,ab - iV^gA^^f ■ (56b) 

Now, the connection on the worldsheet is given by 

TsC = -jl^^ {lBE,C +1EC,B -IBCE ) = j'^^X'',e X^,bC (57) 

and the curvature (either directly or from the Gauss-Codazzi equations) is 

RaBCD = X^j^(jX^^-BD — X'^j^jjX ^_i-BC ■ (58) 

We may therefore, using the identity 

RaBCD = ^-R [lAClBD - lADlBc) (59) 

and the symmetries of N , infer the following useful relations: 

K - \R5f> - -N^i (60a) 
R = -TV^I (60b) 
N^%NES - X^-^^X-^^ {Racbd + X^^^X,,Bc) = \r^ + NigN^E- (60c) 



We can now rewrite the action and its constituent variations as 

1 



/ d'a^ ■>l-e^^R+- 



(^2 + -a,)R' + a^N^ijN^g 



(61) 



and 



9 



SV^ = IV^I^^'S^AB = l^'^X'^ JX^.B (62a) 

5R = -2X«^^<5(X^.^b) + 47Vi'^X^'^^X,,.z3 = -2X^-'^''{5X^),AB - 2RX^^''6X^,d (62b) 

S{NiSN^E) = HSX.y^ABX^^^'^N^E AX^^j^SX^.a^^^nESN^E (62c) 

[using the fact that JF^^X^^ is paraUel to the worldsheet whereas X^'^^ is perpendicular, and (|60|).] 

Noting that / ^—^R oc Xj the Euler characteristic which is a topological invariant, or from the above equations, 
we see that the ai term will not contribute to the equations of motion. The R^ term, on the contrary, gives 

3(7^^X%i?2).^ - 4 {RX^-^^^),j^^ = - AR.aX^-^^^),^ = -AR.abX^''^^ (63) 

where we used the Riemann identity 

{DaDb - DbDa)X^..c = RcdabX''^'^, (64) 



which implies 



at each step. 

The N"^ term gives 



but 



Xt^-^AB ^ = R^^^X^g = (65) 



3{X^'^nESnEE).a + HN^EX'-'^'': 



■AB ' 



atAB \^n:CD vu:B v \A x^u 



CD 



= -\{Ricd),bX''-'^'' (66) 
hence the full equations of motion for the system are: 

^□X'^ = -(2a3 + 4a2)i?;AsX«^^ + ia^^X^^^^ {nESnEE),a + ^c^zN^EX^'^'''' .ab- (67) 
It is possible to express these equations in terms of the extrinsic curvatures. Using ( p4[ ) and 

^^AB - X^B - ^abX% (68a) 

X^^oXf'.AB = 0, (68b) 

it follows that 

Kab = '7^"^.;Ai3 - {i''°X%X';j, - X,,ab 

= -5'=n'tn'<X,^AB = S^'n'tK^AB- (69) 
Contracting the equation of motion with and with X^p gives equations of motion normal and parallel to W: 

^Kfj, = (4a2 + 2a^)R,ABKf'' + ^a^S^'^ Kf K^"" n^^X^,cDAB (70a) 

= Sag (nEEnEE) p + ^asN^El'"''!''^ X^^pX^.,cDAB (70b) 



Note that this latter equation ( |70b| ) is an identity for the unperturbed worldsheet. We have verified that it does 
indeed hold using the lightcone gauge. Now, 



- J""i^„ci3i^f4^.£B - '5™"e™^A;sif„CD + K^CDUJAOJB, (71) 



n'^Xfj^-CDAB = — K^CD-AB — <5™"£imW_Bi4r„cr>;A — ^"^^^ Sirn^^ aK^cD-B 

^^A 

which implies that the equations of motion also read: 

^Kf^ = -(4^2 + 2a^)R,ABKf'' + Aa:iKfKf''{~ K^cd;AB - e,k^BK,cD-A - e^kCOAK.cD-B 

- K^CdK^j^K^EB - £ik(^A;BK^cD + K^cDt^Ai^B) ■ (72) 
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V. CORRECTIONS TO THE MOTION OF TEST TRAJECTORIES 



We now wish to derive the corrections to the motion of three test trajectories: the circular loop, the travelling wave, 
and the helical breather. We have chosen these three examples because they provide three different ways in which to 
observe the rigidity or otherwise of the string. 

The loop trajectory is given by 

X'^ (r, a) = (r, cos r cos tr, cos r sin ct, 0) , (73) 

and collapses to a point after a time period At = 7r/2 = L/A, where L = f da = 2tt is the length of the closed string. 
The extrinsic curvature invariants become singular at this point, hence rigidity would be indicated by a retardation 
of the collapse or a positive correction to the amplitude of the loop. 

A travelling wave on the other hand is a variant of the flat worldsheet where a deformation of arbitrary size and form 
is introduced, the only constraint being that the deformation is a function of only one of the lightcone coordinates, 
a± = cr ± t: 

X^(T,^) = (r,/(T-a),g(r-a),a) (74) 

This has been shown to be a solution of the full field theory hence no correction should be found for this trajectory. 
The helical breather is a time-dependent solution which may be written as 

(r, a) = ^T, v^l — cos r cos cr, \/l — q'^ cos r sin cr, qa^ (75) 

The limit g — > 1 corresponds to the flat worldsheet and the limit q — > to the collapsing circular loop. For intermediate 
q the trajectory is never singular, and the extrinsic curvature peaks at approximately q^^y/l — q^- Rigidity would be 
indicated by a preference for lower extrinsic curvature and hence a negative correction to the amplitude of oscillation. 



A. Corrections to the motion of a collapsing loop 

As a first example of the effects of curvature terms on the motion of freely-moving cosmic strings, we consider the 
collapse of a circular loop. In Cartesian coordinates, the position of this string is given by 

X^(t,o-) = (r, Z(t) cos cr,Z(T) sin cr,0), (76) 

and the normals to the worldsheet are 

= (z, cos cr, sin CT,o) /^l- (77a) 

< = (0,0,0,1). (77b) 

In this case, one obtains uja — 0, and 

K^AB = ( f ) N^-Z^ (78a) 

IC.AB = (78b) 

lAB-(' ) (78c) 



-zz 







Tab = ) (78d) 



^AB = %(\\\- (78e) 



The equation of motion to zeroth order is 



+ 1 - 



^.A = ■■ 7- = 0, (79) 

2^ Z(l-Z2)3/2 ^ ' 
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The general solution to this equation is Z{t) = fccos([r — ro]/fc). Choosing A: = l,ro = as initial conditions one 
obtains the canonical form of the loop trajectory (fz^). Note that for this choice of loop length k = 1 and we can use 
K ov K, interchangeably. 

We now wish to find the corrected solution to order = r^. To do this, we use equation ( [z^ (technically its 
unrescaled counterpart). Since K^ab and uj^ vanish, the RHS of (j?^) considerably simplifies to 

4(a2 + a3)7^;Ai3/C^^ = 32(a2 + ag) secS(T)(7sec^(T) - 6). (80) 

(Note that we have dropped the subscript '2' on the extrinsic curvature.) 

The left-hand side is obtained by varying the trace of K,^ from (|79|), whereby we obtain the equation for 5Z as 



5Z + 2tan(r)5Z - 5Z = 32^ (aa + ag) sec^(T) (7sec2(T) - 6) . (81) 

The solution to this, with initial conditions 5Z{Q) = SZ{Q) = is 

e\ s / 7 1 1 _ 31 , ^ T 



6Z{t) - 32-(a2 + as) (^^ sec5(r) + - sec3(r) + - sec(T) - ^ cos(r) - - sin(T) j . (82) 

Thus, the sign of 0,2 + 03 determines whether the string is rigid or antirigid. This 'loop rigidity parameter' is plotted 
against /3 on figure and its negativity means that the loop is antirigid and tends to collapse faster than in the 
Nambu approximation. This is illustrated on figure H, where we compare the zero-order solution Z{t) — cos(t) with 
the corrected solution Z{t) = cos(t) -|- 5Z{t). 

Note that the approximation breaks down when \K,g\ = 0{rj^)^ i.e. when cos(t) — 0{rs). For example, in the 
illustration of collapse in figure ^, rs — 1/10 is rather large; we would hope that our approximation would be valid until 
the radii of curvature of the worldsheet became close to 1/10, let us say twice the radius of the string: 1/5. Inputting 
\JCg\ = 1/5 gives r ~ 1.1, which does indeed correspond to the point at which the solutions start to significantly differ 
in figure |^. 

B. The motion of a travelling wave 

Now consider a travelling wave, whose position and normals are given by: 

X^(T,a) = (t,/(t- a),.9(T -a), a) (83a) 



<^n^ = (0,5',-/',0)/v//'2+g'2 (83b) 



< ^ = {.r + g'\ /', <?', + .9") l^rTV^ (83c) 



where a prime denotes differentiation with respect to the argument, namely t — a. 
Writing 



one has: 



A(T-^) = /'2+5'2 (84a) 

C(t -cj) = U"g' f'g") /A (84b) 

^2(T-a) = (/V-/V')/VA (84c) 

'^3(r-a) = (/7"+5V')/VA, (84d) 



ujA = Cy\j (85a) 
K^AB = \ ) (85b) 



It is then straightforward to see that all terms on the right-hand side of ( [7^ ) vanish separately, as do the traces of the 
extrinsic curvatures: the travelling wave is an exact solution to (at least) fourth order. 



12 



C. Corrections to the motion of a helical string in breathing mode 



We now consider the string given by the following position functions: 

= (r, Z{t) cos(cr), Z(r) sm{a),qa) (86a) 
= (0, gsin(a), -gcos(a), Z) /^q^ + Z'^ (86b) 

= [z, cos(cr) , sin(cr) , / Vl - Z"^ ■ (86c) 

This string is helical wit h breathi ng q: the limits g — > and 9 — > 1 represent a collapsing loop and a straight string, 
respectively. [Note how (B6a- ^6c ) reduce to (|76|-|77|), and how all pertinent quantities in the equations of motion are 
obtained from these expressions.] 

With this choice of normals, the fundamental forms are: 



lAB 
UJA 



l-Z^ 

-{q^ + Z^) 

qZ /O 1 

1 



Z 
-Z 



Hence the equations of motion to zcroth order become 

Z Z 



(1_ ^2)3/2 (g2_^^2)(l_ ^2)1/2 

This equation admits for general solution 

Z{t) = y/k^ ~ g2 cos (^^^ 

so that, choosing again the initial conditions fc = 1, tq = and calling 

VL{t) =cos2(r) +g2sin2(T), 

we have 



= 0. 



7AB = ^ 



1 

-1 



'1 - g2 
K^AB = q\l ^ sin(r) 



1 

1 



AB 



, . / 1 







n V 1 

The right-hand side of the corrected equation of motion becomes then 



g2 cos(T)f^-9/2 _ p^^-l ^ _ ^_j^-3^ 



where 



/3i = 6 [(a2 + as) + o-iq^] 

P2 = 7{a2 + as) + (38a2 + 22a3)q^ + {7a2 - Sa3)q^ 

= [(7^2 + 5a3) + (7a2 + 2a3)q^] 
/34 = 15g^(2a2 + a3), 



(87a) 
(87b) 
(87c) 
(87d) 



(89) 
(90) 

(91a) 
(91b) 
(91c) 
(91d) 



(92a) 
(92b) 
(92c) 
(92d) 
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and the left-hand side is obtained as before by varying the trace of : 

S {K^a) = ^'^^^SZ + 2n-^^'^{l ~ q^) sin(r) cos(r)5Z + Q-^^^ [q^ - (l - q^) cos'^{t)] 5Z. 
Finally, the corrected equations of motion are: 

6Z + sin(r) cos(r)^Z + l_(l^£h^sZ 



(93) 



-32-yr 



q"^ cos 



(94) 



Although it is possible to find an exact solution to this equation (see appendix for details), it is more instructive to 
consider the quasiflat limit q ^ 1. In this case, equation (P4|) becomes 



5Z + 6Z 



-32 



[(q;2 + as) cos(t) - (2a2 -l- a^) cos(r) sin^(T)] . 



where A is defined by 

q^ = 1- A^. 

The solution 5Z{t) satisfying 5Z{Q) = 5Z{{)) = is then found to be 

e^A^ ' 



5Z 



4 (2a2 -I- 3a3) r sin(r) + (2a2 -I- as) ( 



cos T — COS 



(3t)) 



(95) 
(96) 

(97) 



The corrected trajectory can then be written 
' e4A2 



Z + 5Z 



1 



-(2a2 -I- as) 



4e4A2 , 

1 + (2^2 +3^3) 



1^3 



+ (2a2 + a3)cos(3T) +0(e«). (98) 



The effect of the correction is threefold: First, it alters the frequency of the motion, r [l -I- 4e'*A^(2a2 -f- 3a3)//i] r; 
since (a2 -I- as), 013 < this has the effect of reducing the frequency — a tendency we would be tempted to call 
rigid. Secondly, the amplitude of the oscillation is altered by a factor 1 — e''A^(2a2 -|- a^)/^. This could be either a 
amplification or reduction, depending on whether /3 > 1 or /3 < 1. Finally, a higher frequency oscillation is superposed 
on the motion for /? 7^ 1. If, for simplicity, we take /3 = 1, so 2a2 + = 0, we see that 



Z + SZ = Acqs 



1 



-03 r 



(99) 



i.e. the only effect of the correction is to reduce the frequency of oscillation of the breather, which would seem to be 
unambiguously rigid. 

However, we now observe a curious property: Suppose instead we consider initialising the correction at the instant 
of maximal velocity: 5Z{—Tr/2) = 6Z{—7t/2) = 0, we find 



Z + 5Z 



e^A^ 

1 (2a2 -h9a3) 



sm 



4e4A2 , 
1 + (2a2 -f 3a3) 



Ai 



-(2a2-f a3)sin(3T') + 0(e^) (100) 



(where r' = r + 7r/2). Now note that while the frequency of oscillation is decreased by the same amount, and the 
higher frequency term is still the same, the amplitude is now uniformly increased for all If, as before, we take 
/? = 1, we now find 



Z + 6Z = A 1 



Ai 



as sm 



, Be^A^ 
1 -I as I r 



(101) 



in other words an increase in the amplitude of oscillation accompanies a similar decrease in its frequency. One cannot 
immediately see from this solution whether or not the behaviour is rigid, however, an analysis of the Ricci curvature 
near r' = shows that it is in fact increased — a behaviour consistent with antirigidity! A calculation of the dependence 
of the curvature on general q is presented in the appendix. The algebra is more complicated but the results are the 
same: all helical breathers display both rigid and antirigid characteristics. 

The results of our helical breather calculations therefore appear rather ambiguous. If we wish to identify rigidity 
by the behaviour of a corrected trajectory — whether it increases or decreases curvature — we are forced to calculate 
the effect on the curvature (see appendix) and then the results appear to depend on the initial conditions. What this 
shows is that the 'decrease in curvature' criterion for rigidity is too naive to be reliably applied in all situations. 
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VI. THE QUESTION OF RIGIDITY 



We now want to determine whether a string can be labelled 'rigid' or 'antirigid'. To do this, we use an argument 
based in that of Polyakov , which consists in determining how the action varies under a rescaling of the spacetime 
coordinates X^^ XX ^. Such transformations alter the scale of crinkles of the worldsheet and magnify or reduce 
small-scale structure, hence rigidity would be indicated by a extremum of the energy or the action with respect to 
the rescaling parameter A, as illustrated in figure |ll|. 

Out starting point is the action ( ^o|) , which can be written 




S" = M / d-'ay/^-e^ai / d^ay/^ Mu + e'a2 / dVV^M/, + e^a^ / d^cj^/^ M.^M^j, (102) 



where the matrix My is defined as 



M,j=K^ABKf''. (103) 




As MijMij = Mf^ — 2 det(M), this can be expressed as 

^lA- e^aix + [("2 + as)/! - 2^3/2] , (104) 

with A the area of the worldsheet for the range of {t, ct} being integrated over and x proportional to the Euler 
character of W. Also, 

(105a) 

~ det(Af). (105b) 

We now rescale X^^ AAT'', so that for A > 1 the worldsheet is expanded and for A < 1 it is shrunk. Then, 

lAB -> >?lABj (106a) 
K^AB -> XK^AB, (106b) 

and thus li — s- X~'^Ii, i — 1,2. The shape of the curve S{X) now depends explicitly on the integrals since —S is 
rescaled as 

-S^ X^fiA + e^aix + A^'e^ [(^2 + ^3)/! - 2^3/2] . (107) 

We know that 03, a2 + 03 < 0, and clearly Ii > 0, so in order to determine the shape of ^(A) we only need to 
determine the sign of det(M). For these purposes, we can work in the conformal gauge, ^ab = Vab, and we find 

det(M) = (i^.ooA^ii - i^30oi^,ii)' - 2(/^,ni^3io - K^iiK,_w)^ - 2{K^ooK^io - A^ooif.io)"- (108) 

If we impose the Nambu equations of motion, AT'^ = 0, this determinant is strictly negative. Hence, since ijA is 
positive, we see that S{X) is unbounded, because the the coefficients multiplying A^ and X~^ have opposite signs. 

We must therefore conclude that the string is antirigid. This does not mean that all the trajectories of the string 
exhibit antirigidity, but rather that it is impossible for all trajectories to be rigid. 

Let us illustrate this by considering the trajectories of the previous section. For the collapsing loop, det(M) = 0, 
and as we noticed, 012 + as determines alone the shape of S{X). This simplification came from the fact that the loop 
is flat, and therefore has only one nonvanishing extrinsic curvature. 

For the travelling wave, Mij = 0, which is consistent with our observation of no corrections to this trajectory. 

In the case of the helical breather, both Ma and det(M) are non-zero, so we do not expect results to depend on 
Q!2 + a3- Even though the action is unstable to the scaling of the worldsheet, what is happening with the helical 
breather is that the correction does not always have a non-zero projection on this unstable mode. However, we would 
expect a general correction of the breather to exhibit an instability. 

To sum up: We have reviewed the derivation of the effective action for a U(l) local cosmic string to fourth order in 
the ratio of string width to worldsheet curvature. We presented numerical results calculating the coefficients of these 
fourth order terms. We then derived the equations of motion for the string to fourth order, and calculated corrections 
to a sample of well-known trajectories. We have given a general argument for antirigidity of the cosmic string to 
fourth order, however by reference to our examples have shown that not all trajectories need behave in an 'antirigid' 
fashion — rigidity it appears is rather like a theorem, it may work in special cases, but one needs only find a single 
counterexample to disprove it. 
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APPENDIX: THE HELICAL BREATHER 



Recall from section V C that the helical breather solution to the Nambu action is 

— (r, Acos(T)cos(cr), Acos(T)sin((T),gCT) 

[where A is defined in (|96|)] and that the equation of motion for SZ{t) is (p^: 

A^ ^ ^ , ^- g2 _ ^2 cos^Tt) 
SZ + 2— sin(T) cos{t)SZ + —6Z 

A 

= -32— Acos(T)r2-3(/3^ _ p^n-^ + /33f7-2 _ fj^n-^). 



(Al) 



(A2) 



This can be solved using the method of variation of parameters, giving (after a long and tedious calculation): 



SZ{t) 



16 



3 4 5 6 



)(co*) + AV.„M)-(|L-|L,j|_j|)„,M 



-Aosin(T)tan ^ (gtan(T)) + sin^(r) cos(r) A„r2 



(A3) 



n=l 



Here, ^^(t) = cos^(r) + sin^(r), the /3's were defined in (|9|) and the A's are: 
A 



Al 

A2 

A3 
A4 
As 
Ae 



720/i 



360// 
A 

90^1 
A 

8A 
15^ 



640(3 + l/q^)f3i ~ 360(5 + + l/q^)f32 + 48(35 + 15/q^ + 9/q^ + 5/q^)f33 

-25(63 + 28/g2 + 18/q^ + 12/q^ + 7/q^)f3i 
1920(3 - l/q'^)Pi - 360(15 - A/q^ - 3/q^)P2 + 48(l05 - 25/q^ - U/q^ - 15/g^)/33 

-5(945 - 210/g2 - 136/9"* - HO/q^ - 105/g^)/34 
1920/3i - 360(5 - l/g^)/32 + 48(35 - 6/q^ - 5/q'^)f33 - 5(315 - 49/g^ - 39/9-* - 35/g^)/34 
480^1 - 360/32 + 48(7 - l/q^)f33 - 5(63 - 8/q^ - 7/q^)(3. 
-60^2 + 48^3-5(9- l/(j2) ^4 



6/33 - 5/34 



-T— /J4- 

3fi 



(A4a) 

(A4b) 
(A4c) 
(A4d) 
(A4e) 
(A4f) 
(A4g) 



To use this general solution to investigate the (anti)rigid nature of these helicoidal trajectories, we need to observe 
how the Ricci curvature 



R = K^abK 



AB 



K K = 



2q^Z^ 



2ZZ 



(g2 + Z2)2(l-Z2) (1 - Z2)2(g2 +^2) 



(A5) 
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depends on the correction. For simpHcity, we take = 1, and note that for the background solution 



2A2 



COS' 



(A6) 



Now, suppose we wish to investigate the behaviour of the curvature near a general initial point tq, where SZ{tq) 
SZItq) — 0. Then near tq 



r!8(To) 



f32^ + 



(A7) 



where we have used ( p4[ ) to evaluate SZ{tq), and noted that p4 = for /3 = 1. Now, the combination R5R will 
be neg ativ e if if the magnitude of the curvature is decreased, which corresponds to an intuitive notion of rigidity. 
From (A6) and (A7) we see that this requires: 



[(2 - A^) I]2 - 1] (-A^ + 2A^ - SA-^E^ + ISA^E^ - GA^S"*) > 



(A8) 



where = sin^(ro). Figure |l^ shows the sign of RSR as a function of the two parameters E and A. The shaded zones 
indicate the regions where the string is antirigid. We see that — with the exception of the loop case (A^ — 1) — the 
string admits both rigid and antirigid behaviour for each value of A. 
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FIG. 3. Higi 
for /3 = 1.0. 
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r (units ofr^) 

;s field g (dashed line), polar gauge field q^, (solid line) and radial gauge field qr (in function of r) solving (^) 
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FIG. 4. The parameters ^{l3)/iTrf^ (solid line) and ai{^)/-KTf appearing in the action to fourth order. 
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FIG. 9. The corrected evolution of the quasiflat limit of a helical breather for /? = 1 and e = A = 1/10. The correction 
added to the Nambu solution Z{t) (solid line) is in fact 10* SZ in this figure (dashed line). 
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FIG. 10. The coefficients Q2 + fas (solid line) and 2q2 + as apf^aring in the solution SZ. Note that this last combination 
also appears in the action, and that it vanishes for the critical coupling /3 = 1. 
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FIG. 11. Schematic graph of — ^(A), the variation of the action as the worldsheet is rescaled. The solid line would correspond 
to a rigid string, and the dasfied line to an antirigid string. 




TABLE I. The numerical coefficients appearing in the action to fourth order for some values of the Bogomol'nyi parameter 

(3 IJ'/t^'H^ aij-KT) (xij-K as/ TV 

O 3.272 6.025 77.761 -104.985 
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0.2 2.813 3.347 23.565 -32.835 

0.5 2.314 1.634 4.308 -6.767 

1.0=^ 2.000 0.999 1.069 -2.138 

2.0 1.736 0.642 0.174 -0.722 

10.0 1.278 0.272 -0.097 -0.094 

50.0 0.977 0.141 -0.065 -0.030 

100.0 0.880 0.111 -0.052 -0.023 



''For /J = 1, it can be analytically deduced from the equations of motion that n/irrj^ = iai/irif = 2 and that 2a2 + as = Q. 
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